Toward a structural understanding of turbulent drag reduction: nonlinear coherent 

states in viscoelastic shear flows 
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Nontrivial steady flows have recently been found that capture the main structures of the turbulent 
buffer layer. We study the effects of polymer addition on these "exact coherent states" (ECS) in plane 
Couette flow. Despite the simplicity of the ECS flows, these effects closely mirror those observed 
experimentally: structures shift to larger length scales, wall-normal fluctuations are suppressed 
while streamwise ones are enhanced, and drag is reduced. The mechanism underlying these effects 
is elucidated. These results suggest that the ECS are closely related to buffer layer turbulence. 
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Rheological drag reduction, the suppression by addi- 
tives of skin friction in turbulent flow, has received much 
attention since its discovery in 1947 jl], | S|. For many 
polymer-solvent systems, the pressure drop measured in 
the pipe flow of the solution can be 30 — 50% less than 
for the solvent alone. The central rheological feature 
of drag-reducing additives is their extensional behavior 
in solution: for dilute polymer solutions in particular 
the stresses arising in extensional flow can be orders of 
magnitude larger than those developed in a shear flow. 
This fact is well-recognized; nevertheless the mechanism 
of interaction between polymer stretching and turbu- 
lent structure is not well-understood and the goal of the 
present work is to attempt to shed light on this interac- 
tion. 

A key structural observation from experiments and di- 
rect numerical simulations (DNS) of drag-reducing so- 
lutions is the modification of the buffer region near the 
wall f§ || g, |, | 0. It has long been known that the 
flow in this region is very structured, containing stream- 
wise vortices that lead to streaks in the streamwise ve- 
locity JO] ]; these structures are thickened in both the 
wall-normal and spanwise directions during flow of drag 
reducing solutions 0, M. Because of its importance in 
the production and dissipation of turbulent energy p| , 
any effort to mechanistically understand rheological drag 
reduction should address this region. 

To better understand the effect of the polymer on 
the buffer layer, we wish to study a model flow that 
has structures similar to those seen in this region but 
without the full complexities of time-dependent turbulent 
flows. Fortunately, a family of such flows exists, in the 
recently-discovered "exact coherent states" (ECS) found 
by computational bifurcation analysis in plane Couette 
and plane Poiseuille flows ||, [0| g |l5|, ^BJ. These are 
three-dimensional, traveling wave flows (hence steady in 
a traveling reference frame) that appear via saddle-node 
bifurcations |35} | at a Reynolds number somewhat below 
the transition value seen in experiments |17L 18 . The 



structure of the ECS captures the counter-rotating stag- 
gered streamwise vortices that dominate the structure in 
the buffer region. From the dynamical point of view, 
there is evidence that these states form a part of the 
dynamical skeleton of the turbulent flow: i.e., they are 
saddle points that underlie the strange attractor of turbu- 
lence [|l9| . Finally, the nonlinear self-sustaining mech- 
anism underlying these states has been elucidated J2lj . A 
perturbation of the base flow in the form of streamwise 
vortices redistributes the streamwise momentum of the 
flow. This redistribution creates spanwise fluctuations 
in the streamwise velocity, the "streaks". The spanwise 
inflections in the streamwise velocity profile lead to a 
three-dimensional instability that develops into staggered 
nearly-streamwise vortices that regenerate the streaks. 
Because the ECS capture the structures of the buffer 
region and are mechanistically well-understood, we be- 
lieve that they provide an excellent simplified, yet still 
exact, model flow for studying polymer drag reduction. 
The leading order effect of viscoelasticity on the ECS is 
therefore the focus of the present study. 

To begin, we briefly describe a general result relating 
polymer stretch to flow kinematics. For a trajectory in 
a flow field, the Liapunov exponents give the Lagrangian 
time-averaged rate of stretch of material lines. If the 
largest Liapunov exponent, cr max , is positive, the flow is 
extensional on average. In particular, for homogeneous 
turbulence, the expected value of a max is positive |||] 
and we show below that this is also the case for the ECS. 
Now consider the dynamics of a Hookean dumbbell model 
of a polymer in a flow field. The end-to-end vector q of 
the dumbbell evolves in the flow field, v, as: 
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where D/Dt is the time derivative evaluated on a fluid 
element, A is the stress relaxation time for the dumbbell 
and £(i) is the random Brownian force. Noting that an 
infinitesimal material line satisfies the same expression 
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but with A -1 and £ set to zero, it is straightforward to 
show that Hookcan dumbbells will stretch indefinitely in 
a flow if and only if Xa max = We a > s, where We a 
is a Weissenberg number based on a max . This result 
is a specific statement of an idea that originated with 
Lumley |23j (see also g |§ |7j)- The computations 
below show the importance of We a in determining the 
effect of polymers on coherent structures. 

We study here the effect of polymer on the exact co- 
herent states that arise in a variant of plane Couette flow 
|[4| . Denoting the streamwise direction as x, the wall- 
normal direction as y, and the spanwise, or vorticity, di- 
rection as z, we consider a flow with boundary conditions 
^ = l t v y = ^ = at y = ±1. The characteristic 
velocity, U , and the half-height of the channel, /, have 
been used to scale the velocity and positions, respec- 
tively. These "constant vorticity" boundary conditions 
provide an advantage over no-slip conditions in that they 
allow us to model only the buffer region in our domain by 
eliminating the viscous sublayer. (Exact coherent states 
found using no-slip BCs jl6) show a qualitatively iden- 
tical vortical structure, only offset from the wall by a 
small region comprising the viscous sublayer.) Periodic 
boundary conditions are applied in the streamwise and 
spanwise directions. For this study, the wavelength of 
the structures in the streamwise and spanwise directions 
are fixed at £ x = 27r/0.40 and £ z = 27r/1.0, respectively 
(or 165 and 66, if expressed in wall units at a Reynolds 
number of 110). For this flow, a trivial (Couette) base 
state exists, v x {y) = y; the maximum mean velocity for 
the ECS is significantly reduced compared to the base 
state velocity due to the enhanced transport of momen- 
tum @. 

In our formulation, time, t, is scaled with l/U, and 
pressure, p, with pU 2 , where p is the fluid density. The 
stress due to the polymer, x p , is nondimensionalized with 
the polymer elastic modulus, G = i] p /X, where rj p is the 
polymer contribution to the zero-shear viscosity and A 
is the relaxation time for the polymer. The momentum 
and mass balances are 
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where r] s is the solvent viscosity, El = A(rj s + rj p )/pl 2 and 
(3 = T) s / (r] s + rip). The Reynolds number, Re, is based on 
the total viscosity, Re = pUl/(r) s +r) p ). 

We calculate the polymer stress with the commonly 
used FENE-P model p8[ |, which idealizes the polymer 
molecules as bead-spring dumbbells with finitely extensi- 
ble springs. With this model, the non-dimensional struc- 
ture tensor a (= (qq), where () denotes ensemble aver- 
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where We - 



is the Weissenberg number based on the 



wall shear rate and 6 is proportional to the maximum ex- 
tension of the dumbbell: trot cannot exceed 6. A simple 
measure of the importance of extensional polymer stress 
is the magnitude of the parameter Ex = | . In uniax- 
ial extension with extension rate i, Ex = 1 implies that 
Tp = t v as e — ► oo where t v is the solvent contribution 
to the stress. The polymer can significantly affect the 
flow only when Ex > 1 . In the present flow this parame- 
ter also represents the maximum ratio of polymer stress 
to Reynolds shear stress (the flux of i-momentum due 
to fluctuations in y- velocity) ; the scaling theory of Tabor 
and de Genncs 
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51j treats the regime Ex ^> 1. 
The governing equations are solved through a Picard 
iteration. A given velocity field is first used to calculate 
the polymer stress tensor, t p , by time-integrating Eq. ^ 
until a steady state is attained. For the new r p , a steady 
state of the momentum and continuity equations is found 
by Newton iteration. The resulting velocity field, is used 
to compute the new t p , and the process is repeated un- 
til the velocity field converges. Equations |||| are dis- 
cretized as in [Q, using a Fourier-Galerkin formulation 
with typically a 7x19x7 grid. Equation || is discretized 
with a Fourier-pseudospectral method, typically with a 
32x32x32 grid, and time-integration performed with an 
Adams-Bashforth method. To achieve numerical stabil- 
ity, a small diffusive term is added to Eq. || (cf. ||) and 
integrated with a Crank-Nicholson scheme. 

Before presenting the effects of the polymer on the 
ECS, we recall the result that We a > 1/2 implies large 
stretch of polymer chains. For the Newtonian ECS at 
Re — 110 on the lower branch of the bifurcation diagram 
(see Fig. [l]), the velocity field is very nearly ergodic, with 
Umax ~ 0.030. The condition We a > 1/2 thus trans- 
lates into We > 17 for large polymer stretch, and for 
Ex = 0(1), will define the onset condition for the poly- 
mer to begin to strongly affect the flow field. In DNS 
of a FENE-P fluid in plane channel flow, Surcshkumar, 
et al. J9|, found no drag reduction at We = 12.5 and 
significant drag reduction at We — 25; another recent 
DNS study §§] places the onset value at We w 20. This 
close correspondence between the onset condition pre- 
dicted from the ECS kinematics and that found by DNS 
strongly suggests that the ECS model captures the es- 
sential structure of the buffer layer. 

Fig. shows how the addition of polymer stress af- 
fects the bifurcation diagram (i. e. the locus of steady 
state flows) for El = 0.20 {We a « 2/3) and Ex = 1 
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FIG. 1: Newtonian and viscoelastic (El = 0.20, (3 = 0.97) 
bifurcation diagrams. Each curve represents a locus of steady 
state ECS flows; the leftmost point on each curve is the posi- 
tion of the saddle-node bifurcation. 



to 3. On the j/-axis of the diagram is the maximum 
value of the streamwise- and spanwise-averaged stream- 
wise velocity {(v x )). (The trivial Couette solution is at 
max((v x )) = 1.) When Re attains a certain value that 
we denote Re sn , two new steady solutions appear in a 
saddle-node bifurcation. For small values of Ex, Re sn 
decreases compared to the Newtonian value, but once the 
polymer stress begins to exceed the viscous (Ex ms 1.5), 
Re sn increases above the Newtonian value - the pres- 
ence of the polymer is suppressing the ECS. Qualitatively 
identical behavior is observed experimentally in the onset 
Reynolds number for turbulent pipe flow 0. Note that, 
for a given velocity U and total viscosity 77, the increase 
in critical Reynolds number corresponds to an increase in 
the characteristic length scale of the coherent structure, 
again consistent with experimental observations ||, ||. 
Finally, we see that ((v x )) becomes larger for the vis- 
coelastic flows than for the Newtonian - drag reduction 
occurs. 

To examine more closely the effect of the polymer 
stress on the velocity fields, Fig. ^ shows results at con- 
stant Re while varying El, or, equivalently, We a (based 
on <r max = 0.030). Here we plot ma,x({v'y)), where the 
prime denotes the fluctuating part of a quantity, here 
wall-normal velocity. At Ex — 1, after an initial in- 
crease, max({v'y)) decreases below the Newtonian value 
and eventually saturates, as the polymer stress asymp- 
totes at high We G to a fixed value relative to the viscous 
stress. In this case, the polymer stretch becomes nearly 
uniformly large throughout the domain. The decrease in 
wall-normal velocity with We a is even more drastic as 
the extensibility parameter Ex increases. Similar trends 
are seen in the streamwise enstrophy and Reynolds shear 
stress. Since Ex is related to the extensional viscosity 



of the viscoelastic solution, these results show the im- 
portance of extensional stresses in affecting the ECS. 
The spatial maximum of tra, which is proportional to 
the square of the polymer extension, is also presented 
in Fig. |^. Finally, note that the majority of both the 
polymer stretch and the change in wall-normal velocity 
occurs in the range 0.1 < We a < 1.0. In contrast to 
the decreases in wall normal fluctuations, streamwise en- 
strophy and Reynolds shear stress for We a > 0.1, the 
streamwise fluctuations ((v x )} are found to increase. All 
of these trends are observed in DNS and experiments 
ill- 
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FIG. 2: Wall normal velocity and maximum polymer stretch 
vs. El and We a , lower branch solutions, Re = 110, P — 0.97. 

Turning now to how, and why, the flow structure is 
changed by the polymer dynamics, figures [|a-b show (a) 
the streamwise velocity v x at y — 0.875 (where the max- 
imum in the polymer stress occurs) for the Newtonian 
lower branch solution at Re =110 and (b) the difference 
v x ,ve — Vx.n between the viscoelastic (VE) and New- 
tonian (N) solutions. Here we see the "streak" (white 
ribbon) and - by adding the pattern of fig. ||b to that 
of ||a - that this streak is "straightened out" by the vis- 
coelasticity. Figure ||c shows a contour plot of f x , the 
x-component of the force / exerted by the polymer on 
the fluid, (/ = (1 - f3)Re' 2 El~ 1 V ■ t p ), at y = 0.875. 
This force is significantly negative and corresponds spa- 
tially to the region where fluid elements are leaving the 
streak to move into the vorticity-dominatcd regions. The 
y and z components of the force have also been examined; 
they are smaller than the x-component but are clearly 
seen to work directly against the vortex motions: e.g. 
where v y is highly positive in an upwelling, f y is highly 
negative. This behavior is also seen in the buffer layer 
structures in the DNS study of a drag-reducing polymer 
solution by De Angelis et al. The origin for this 

structure of the polymer force field becomes apparent on 
examination of polymer stresses along fluid trajectories: 
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FIG. 3: Streamwise velocity for the Newtonian (N) and vis- 
coelastic (VE) solutions and polymer force at y = 0.875. 
Lower branch solutions at Re = 110 and for the viscoelastic 
solution El — 0.15, Ex = 3, and /3 = 0.97. (a) v Xi n (range: 
0.0 (black) — 0.889 (white)) (b) v x ,ve - v x , N (range: -0.0395 
(black) — 0.0395 (white)) (c) f x (range: -0.00413 (black) - 
0.00413 (white)). By symmetry, identical behavior with op- 
posite signed velocity is occurring at y — —0.875. 



polymer molecules stretch in (or moving into) the streak 
regions, remaining highly stretched until they begin to 
leave the streak. As molecules move from the streak into 
and around the vortices, they relax. The spatial gradi- 
ents in stress accompanying this relaxation work against 
the vortices, "unwinding" them. This vortex suppression 
leads to collapse of the self-sustained process - or more 
precisely to a shift of the process to larger scales - and 
thus to drag reduction. 

To summarize, we list several points of agreement be- 
tween our results and observations from DNS and exper- 
iments in fully turbulent flow, namely: (1) the ECS bear 
a strong similarity to the structures observed or educed 
from structural studies of the buffer layer and apparently 
underlie its dynamics; (2) for O(l) values of Ex, the onset 
Weissenberg number for drag reduction predicted from 
the ECS kinematics agrees closely with DNS results; (3) 
the effects of viscoelasticity on the Couette flow ECS are 
very similar to those observed in the buffer layer: (a) wall 
normal velocity fluctuations are suppressed and stream- 
wise ones enhanced, (b) Reynolds shear stress decreases, 
(c) streamwise vorticity decreases, (d) the velocity fluctu- 
ations and polymer force are anticorrelated, and (e) drag 



is reduced. Finally, at fixed U and r\ the upward shift in 
the onset Reynolds number corresponds to an increase in 
length scale for the structures, again consistent with ex- 
periments. These successes show that studying the ECS 
holds promise for capturing the essential physics of drag 
reduction. Indirectly, they also validate the view that the 
ECS underlie Newtonian turbulence, because the effects 
of polymers on the ECS so closely mirror their effects on 
full turbulence. 
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